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ABSTRACT 



Spectral lags (xiag) are deduced for 1437 long (T90 > 2 s) BATSE gamma-ray bursts (GRBs) 
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with peak flux Fp > 0.25 photons cm s , near to the BATSE trigger threshold. The lags are 
modeled to approximate the observed distribution in the Fp-xiag plane, realizing a noise-free 
representation. Assuming a two-branch lag-luminosity relationship, the lags are self-consistently 
corrected for cosmological effects to yield distributions in luminosity, distance, and redshift. The 
results have several consequences for GRB populations and for unified gamma-ray/afterglow 
scenarios which would account for afterglow break times and gamma-ray spectral evolution in 
terms of jet opening angle, viewing angle, or a profiled jet with variable Lorentz factor: 

A component of the burst sample is identified - those with few, wide pulses, lags of a few 
tenths to several seconds, and soft spectra - whose Log[ISr]-Log[Fp] distribution approximates a 
-3/2 power-law, suggesting homogeneity and thus relatively nearby sources. The proportion of 
these long-lag bursts increases from negligible among bright BATSE bursts to ~ 50% at trigger 
threshold. Bursts with very long lags, ~ 1-2 < xiag (s) < 10, show a tendency to concenfrate near 
the Supergalactic Plane with a quadrupole moment of ~ -0.10 D 0.04. GRB 980425 (SN 
1998bw) is a member of this subsample of ~ 90 bursts with estimated distances < 100 Mpc. The 
frequency of the observed ultra-low luminosity bursts is ~ _ that of SNe Ib/c within the same 
volume. If truly nearby, the core-collapse events associated with these GRBs might produce 
gravitational radiation detectable by LIGO-II. Such nearby bursts might also help explain 
flattening of the cosmic ray spectrum at ultra-high energies, as observed by AGASA. 

In a regime limited by BATSE sensitivity, 10~^'^< L53 < 10"^'^, the model lags predict a 
power-law scaling relation for the ultra-low luminosity GRBs, dNgen/dL ~ flatter than 
expected (dNsen/dL~L^^^) if viewing angle with respect to the jet axis alone governed perceived 
luminosity. For high-luminosity bursts, in the volume-limited regime z < 2 and 10"^'^ < L53 < 
10^'^, BATSE samples through the distribution and dNyoi/dL ~ L~^'^, similar to expectations for 
viewing angle scenarios (dNyoi/dL ~ L ). However, in the latter case if the luminosity decreases 
off axis, L ~ Oyiew^ > 0), then overproduction of low-luminosity bursts cannot be avoided. 
Thus, a completely relativistic kinematic explanation for the djaiamic range in GRB luminosities 
is not favored. The variable beaming fraction scenario, with constant luminosity across the jet 
cone, can fit the high-luminosity bursts with a fairly flat distribution in jet-cone solid angle, 
dN(Qjet)/dQjet D ^61*^^5 for the ultra-low luminosity bursts a distribution that increases is 
required, dN(Qjet)/dQjet D ^jet^''^- Jets with variable luminosity proflles viewed at a range of 
angles can also reproduce the observed luminosity distributions, such that L ~ Oyiew^'^ and L ~ 
6view^'^5 for high and ultra-low luminosity regimes, respectively. For both the beaming fraction 
and profiled jet scenarios, a large fraction of the SN Ib/c population in the Universe would be 
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required to produce the GRBs at cosmological distances, whose rate is estimated to be 1-few 
xlO^ yr~^. The modeled redshift distribution for GRBs peaks at z ~ 10, with large uncertainty. 

Subject headings: gamma rays: bursts, afterglows, beaming — supemovae: type Ib/c — general 
relativity: gravitational waves — cosmic rays: ultra-high energy 
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1. INTRODUCTION 



Recent theoretical and observational works have appeared which reveal a new understanding 
of the y-ray burst (GRB) energy paradigm compared to the original picture developed over the 
preceding 4—5 years. The new studies conclude that the actual dynamic ranges of energy release 
and luminosity are relatively narrow, rather than spanning factors of a few hundred as implied if 
GRB emissions were isotropic. (Much larger inferred djaiamic ranges [~ 10^] apply if bursts like 
GRB 980425 should be considered within the same framework.) All the explanations involve 
anisofropic ejecta. However, they invoke three distinct physical mechanisms, or combinations 
thereof, and marshal the observational facts to support the differing interpretations. A perceived 
range in luminosity, and therefore energetics, could result from: variation in jet cone opening 
angle while maintaining constant total energy release; variation in viewing angle alone; or a 
profiled jet where the Lorentz factor (F) decreases as viewed off the jet axis. In each scenario, 
evidence has been cited that unifies the y-ray spectral-temporal behavior with the afterglow 
temporal behavior through either pure relativistic kinematics or the jet dynamics. Since each 
explanation may be viewed as economical - each realizes reduced and narrow ranges in y-ray 
energy release and luminosity - observations may provide the key to distinguishing between the 
possible scenarios. There are > 2700 BATSE triggered bursts but only two dozen GRBs with 
associated redshifts, and so y-ray observations may help distinguish between the predictions of 
the three scenarios vis-a-vis the GRB luminosity disfribution as inferred from spectral lags. 

The djoiamic range in energy for cosmologically distant GRBs is ~ 3 x 10^^ - 10^"* ergs, 
assuming isotropic emission. By analysis of break times in the decays of GRB afterglows. Frail 
et al. (2001) have inferred a disfribution of initial jet opening angles, inversely correlated with 
total y-ray energy and luminosity, with most bursts' ejecta exploding into narrow cones, where 
the smallest half angle 9jet > 2-30 The analysis embodies a crucial assumption - that luminosity 
(i.e. F) is constant across the jet's cone. The derived beaming fraction, fb D S2jet/4jt « _ 9jet , 
exhibits a djoiamic range commensurate with observed y-ray luminosities and total energies. The 
true y-ray energy release distribution is then inferred by Frail et al. to be much narrower, less than 
a decade, and centered near ~ 5 x 10^'^ ergs. The temporal signature of afterglow evolution was 
enunciated by Rhoads (1997; 1999) for constant F within the jet cone. When the ejecta's 
Lorentz factor decreases below ~ 6jet ^ manifesting a break in the power-law decay, at that time 
spherical and jet geometries become distinguishable; frans verse expansion of the jet adds to the 
steepening decay. Panaitescu & Kumar (2001: PK) perform an analysis similar to Frail et al. 
The three afterglows in common to the two freatments - for which there are also BATSE data - 
are on the high end of the GRB luminosity disfribution, and the derived sets of opening angles do 
not differ markedly. Figure 1 illusfrates derived beaming fraction and specfral lag for eight bursts 
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with associated redshifts. The two quantities appear correlated over a dynamic range of ~ 50 in 
both coordinates, as would be expected since spectral lag and beaming fraction are inversely 
correlated with observed luminosity. So even though the distance scales from the source (Piran 

13-14 

2001) for the two phases are disparate in the intemal+extemal shocks paradigm (~ 10 cm 
compared to > 10 cm), the correlation evident in Figure 1 implies that the y-radiation djoiamics 
are related to jet dynamics. An adequate parametric fit of the form 

Tiag = -0.127+ 1.11 fb or 
fb = {['Ciag + 0.127]/ 1.11}^°^^ fb^O.05 (la) 

was found for the six bursts with 0.835 s z s 1.619, using the geometric means of beaming 
fraction for the three bursts in common to Frail et al. and PK. GRBs 000131 and 971214 were 
not included in the fit since redshift corrections to their spectral lag measurements could prove to 
be substantial (Morales et al. 2002). The redshift correction factors applied to the lags are small 
(the correction procedure is discussed in section 2.4). The Xiag- fb fit is indicated by the solid 
line. The fit nearly asymptotes near fb « 7.5 x 10""*, equivalent to 6jet « 2.3Q necessarily 
consistent with the minimum jet cone radii of 2-30 which PK and Frail et al. discuss. 

The distribution of BATSE GRB lags described in the next section extends up to 10 s. This 
is a factor of 50 longer than the longest lag (GRB 970508: 0.2 s, uncorrected for cosmology) 
plotted in Figure 1 . Thus, if there were one continuous relationship between lag and beaming 
fraction for GRBs, one would need to construct a continuation of Eq. (la) that asymptotes to fb ~ 
1 at lags _ 10 s. The toy model plotted with a dashed line in Figure 1 achieves the continuity and 
limit requirements: 

"flag = - 0.84 + exp(3. fb) , or 
fb = ln[Tiag + 0.84] 13. fb D 0.05 . (lb) 

Such a change in fijnctional form might reflect a modification of the single power-law form of the 
lag-luminosity relationship in Norris, Marani, & Bonnell (2000: NMB), i.e., to accommodate the 
ultra-low luminosity GRB 980425 - currently an example of one. Salmonson (2001) argues for a 
separate lag-luminosity branch for bursts with low Lorentz factors. However, the proposals by 
Nakamura (1999) and Salmonson (2001) for GRB 980425 involve a large off-jet-axis viewing 
angle, rather than a large jet opening angle. 

Several other investigators have recently discussed alternative interpretations of the afterglow 
breaks that do not necessarily require a range in 9jet. One alternative to the variable beaming 
fraction scenario, invoking viewing angle, has been advanced by several authors, mostly in the 
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context of the low luminosity of GRB 980425 (Wang & Wheeler 1998; Hoflich, Wheeler, & 
Wang 1999; Nakamura 1999; Woosley, Eastman, & Schmidt 1999; Sahnonson 2001) and 
extended to explain the lag-luminosity relationship in general (Salmonson 2000; loka & Nakamura 
2001; Salmonson & Galama 2001). In essence, softer-spectrum, lower-luminosity, longer-lag 
bursts will be observed as the observer's viewing angle off the jet axis increases, merely due to 
relativistic kinematics. 

The third possibility - a profiled jet with T decreasing off axis - is also naturally expected 
fi-om simulations (see e.g., MacFadyen, Woosley, & Heger 2001). Salmonson (2000) and 
Salmonson & Galama (2001) discuss this scenario as well, and argue that the observer's viewing 
angle could explain the apparent djoiamic range in luminosity, pulse spectral evolution, and 
afterglow break times. The profiled jet model is developed in more detail in Rossi, Lazziti, & 
Rees (2002), and discussed by Zhang & Meszaros (2002) in terms of predicted dependence of 
luminosity on angle with respect to the jet axis. 

Independent evidence supports the picture of a small dynamic range in GRB energies. Piran 
et al. (2001) have shown from analyses of X-ray afterglows that inferred total release energies 
span less than a order of magnitude - irrespective of assumptions concerning magnetic field 
strength, electron energy distribution, and external medium. Hence, it is necessary to consider 
carefiiUy the three explanations for a small range in y-ray energy release and luminosity. 

So far, theoretical interpretations of burst pulses reproduce many, but not all of the observed 
pulse behavior. The empirical schema to be elucidated is that at y-ray energies, bursts consist of 
pulses organized in time and energy (Norris et al. 1996). The rise-to-decay ratio is unity or less; 
as this ratio decreases, pulses tend to be wider, the pulse cenfroid is shifted to later times at 
lower energies, and pulses tend to be spectrally softer. However, since ~ 80% of pulses overlap 
in bursts, an extensive parameter - such as spectral lag - must serve as a surrogate measurement 
for the spectral dependence of average pulse shape within a given burst. Variability parameters 
correlated with GRB luminosity have also been reported (Fenimore & Ramirez-Ruiz 1999; 
Reichart et al. 2000). (Specfral lag, variability, and pulse shape may evolve during a burst, but 
the degree of evolution is yet to be well quantified for determination of trends in a large sample.) 
The canonical theoretical schema (e.g., Piran 2001) is that these pulses of y-ray emission are 
produced by the "internal" shocks of colliding pairs of relativistic shells ejected by the central 
engine. The major timescales are accounted for in this picture: overall burst durations and 
intervals between pulses are related to the activity of the central engine, and pulse widths are 
related to the shell collision timescale (Rees & Meszaros 1994; Ramirez-Ruiz & Merloni 2001). 
Details of pulse shape theory are incomplete. A shift in pulse peak as a function of energy is not 
produced by pure relativistic kinematics in the standard picture where the observer is within the 
jet's half angle, 6jet (Fenimore, Madras, & Nayakshin 1996; see also, modeling by Panaitescu & 
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Meszaros 1998). Outside of this cone, the observed pulse is broadened as viewing angle 
increases, and peak shift may be effected by assuming an inhomogeneous Lorentz factor over the 
face of the jet, but this "beyond the jet cone" emission quickly falls below the y-ray portion of 
the spectrum (see Nakamura 1999; loka & Nakamura 2001; Sahnonson 2001; Salmonson & 
Galama 2001; and Kumar & Panaitescu 2000 for off-axis afterglows). Also, the observed locus in 
time of pulse intensity versus peak in v[]F(v) - essentially defining pulse shape as function of 
energy and time - is not reproducible by relativistic kinematics, even when combined with 
variations in T or v[]F(v) across the blastwave emitting surface and with variable shell thicknesses; 
the implication is that pulse spectral evolution requires attention to in situ cooling details 
(Soderberg & Fenimore 2001). And, while spherical blastfront curvature relates the Lorentz 
factor to the observable emission cone (6^ ~ l/F) by the relativistic trigonometry - jdelding a 
simplistic, energy-independent form for pulse width, Atpuise ~ Rshock/2cr^ - only recent 
discussions have connected jet djoiamics to luminosity, pulse spectral lag and/or variability 
(Zhang & Woosley 2002; Rossi et al. 2002). 

It has been suspected that a subclass of ultra-low luminosity, very long-lag, soft- spectrum, 
nearby GRBs exists - GRB 980425 being the well-known exemplar - that could be produced by 
a version of the coUapsar model (MacFadyen, Woosley, & Heger 2001). Beyond the deduction 
that such bursts should have low Lorentz factors (F ~ few: Kulkami et al. 1998; Woosley & 
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MacFadyen 1999; Salmonson 2001) compared to F ~ 10 -10 for the high-luminosity bursts at 
cosmological distances, additional explanations advanced for low observed luminosity include 
viewing angle, profiled jets, or much wider jet opening angles. Kulkami et al. and Wieringa, 
Kulkami, & Frail (1999) infer the latter for GRB 980425, based on the conclusion that the radio 
emission was not strongly beamed. Thus GRBs with very disparate observed luminosities, 
ranging over ~ 10^^-10^^ ergs s"\ may manifest a range of jet opening angles, or they may 
produce profiled jets with variable F which are viewed from different angles. The whole GRB 
population might be unified in one of these senses. 

Also, it is clear that pulse frequency of occurrence within a given burst depends on peak flux. 
Norris, Scargle, & Bonnell (2001) showed in a brightness-independent analysis (to peak flux Fp > 
1.3 photons cm"^ s"') that BATSE bursts with relatively long lags tend to have just a few 
significant, wide pulses and that such bursts are observed preferentially at lower peak fluxes. 
This result was presaged by the burst "complexity parameter" of Stem, Poutanen, & Svensson 
(1999), whose analysis of average GRB profiles as a fimction of peak flux suggested an admixture 
of a larger fraction of simple bursts near the BATSE trigger threshold, Fp ~ 0.25 photon cm"^ s"^ 
A hint that these simple, dim bursts come from sources at low redshifts is provided by the only 
exemplar with known distance, GRB 980425 / SN 1998bw. Its ultra-low luminosity may be 
attributable to membership on a steeper branch in the lag-luminosity "HR diagram" of GRBs (see 
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Salmonson 2001; also Kulkami et al. 1998). The steeper slope of this second branch may be 
related to mildly relativistic outflow; if the initial ejecta have I^' > 6jet, then different behavior is 
expected for the y-ray and immediate afterglow phases. Central questions conceming these 
bursts with relatively simple temporal profiles include their observed frequency, typical beaming 
fraction and luminosity, and the implied volume sampled by BATSE. 

The program here is to estimate luminosities and distances from measured specfral lags for 
long (T90 > 2 s) BATSE bursts, and derive constraints on the GRB population and jet 
mechanism. In Norris, Marani, & Bonnell (2000: NMB) an anti-correlation between y-ray 
specfral lag and GRB peak luminosity was reported, based on six bursts with associated 
redshifts. The relationship is roughly determined as L53 Xiag^'^^. The frend is qualitatively 
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sfrengthened by the addition of two bursts: GRB 991216, Xiag ~ 10 ms, Lpeak ~ 6 x 10 ergs s ; 
and GRB 000131, Xiag ~ 5 ms, Lpeak ~ 10^^ ergs s"\ The original lag values (and the two listed 
here) reported for this relationship were corrected for time dilation, but not for spectral redshift. 
Direct methods for specfral correction involve attempts to deredshift the specfra and measure 
lags, or to construct interpolation tables with lag measured between several pairs of energies 
(Morales, Norris, & Bonnell 2000). Here, only bursts within a range of (1+Zinax)/(1+Ziiiin) =1.4 
are used to define a relationship between lag and beaming fraction, since shifts of the specfral 
energy disfribution for GRB 971214 (z = 3.14) and GRB 000131 (z = 4.5) are large compared to 
the median redshift (z ~ 1) for the balance of the sample. An indirect method for approximating 
specfral redshift correction used the average width of the auto-correlation fijnction versus energy, 
W<ACF> D E"''^'^ (Fenimore et al. 1995). Similarly, the average pulse width as a function of 
energy can be utihzed, W<puise> D E ^-^^ (Norris et al. 1996), as is done here. Both fimctions were 
derived using temporal profiles of bright bursts and, strictly speaking, are applicable over a 
relatively narrow range in redshift. Thus, calibration of the lag-luminosity relation and necessary 
specfral corrections will benefit from denser observations in frequency, over a broader dynamic 
range, and from a larger burst sample with redshifts. 

Section 2 describes the extension of specfral lag analysis for BATSE bursts to near the 
instrument's trigger threshold. The lags are modeled, including noise, to reproduce the observed 
scatter plot in the Xiag-Fp plane. The Log[N]-Log[Fp] disfribution and sky disfribution of the 
long-lag burst subsample are examined. Then assuming a two-branch lag-luminosity relation, 
approximate correction is made for cosmological effects, and the model lags are used to estimate 
disfributions of luminosity, redshift, distance, and beaming fraction. In section 3 the modeled 
luminosity disfribution is compared with expectations for the three jet scenarios described above. 
In section 4 the results are summarized. Some topics in high-energy asfrophysics related to GRB 
studies are briefly discussed, including: the GRB-SN Ib/c connection; rates for observable 
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gravitational radiation from GRB within the coUapsar scenario; the Swift GRB yield; and finally 
the possibility of ultra-high energy cosmic rays from nearby GRBs. 

2. SPECTRAL LAG ANALYSIS 

The first step is to measure specfral lags for nearly the complete sample of BATSE triggered 
GRBs with T90 durations > 2 s. Then using a small number of parameters, the lag-peak flux 
distribution is modeled, in the process realizing a noise-free representation for the lags. An 
iterative procedure is used to correct these modeled lags for the extrinsic cosmological effects, 
while simultaneously unfolding approximations for the GRB luminosity and redshift 
disfributions. This procedure requires an assumption about the correction required for redshift of 
burst specfra, and so yields parameter-dependent distributions, rather than unique results. 

Beyond selecting bursts with T90 > 2 s, the sample is fijrther restricted by requiring Fp 
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(50-300 keV) > 0.25 photons cm s (measured on 256-ms timescale), and peak intensity (PI) > 
1000 count s'^ (> 25 keV). Background fits and burst regions were defined, and peak fluxes and 
durations measured following the same procedures described in Norris et al. (1996) and Bonnell 
et al. (1997). Starting with an available sample of 2699 BATSE bursts, 2024 survived with 
usable, concatenated DISCLA, PREB, and DISCSC data, and with satisfactory background fits; 
1474 of these were measured to have T90 > 2 s. Twenty bursts had 700 < PI (count s"^) < 1000; 
within this lowest count rate range the lag analysis becomes increasingly less useful. 

2. 1 Cross-Correlation Lag Analysis 
A cross-correlation analysis of BATSE channels 1 (25-50 keV) and 3 (100-300 keV) was 
performed, as described in detail in NMB. The peak in the cross-correlation function (CCF) was 
taken as the measure of specfral lag. Only three important modifications to the original 
procedure were implemented. First, the native 64 ms data were binned to 128 ms (256 ms) 
resolution for bursts with PI below 7000 (1400), to facilitate location of the cenfral peak in the 
CCF at low intensities. Second, the 101 realizations per burst with added Poisson-disfributed 
noise (per energy channel), used to estimate statistical errors for the CCF, were restarted with a 
longer fitted range near the peak of the CCF if the fit was concave up. Third, a frick was used to 
eliminate a problem in the IDL poljoiomial fimction which can occur when the independent 
variable range is too narrow compared to that of the dependent variable: the time coordinate was 
expanded by a factor of ten, and the fitted polynomial coefficients were adjusted accordingly. 
The CCF was computed over the portion of a burst temporal profile extending out to the furthest 
points attaining half of the peak intensity. Utilizing outlying, lower intensity portions of a burst 
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results in larger lag errors for the majority of bursts treated here. A cubic fit to the CCF was 
employed to accommodate the asjmimetric nature of bursts on all timescales (Nemiroff et al. 
1994). Seven bursts were eliminated when, during the lag measurement process, the data were 
found to be corrupted by either very intense intervals (and therefore counter overflow) or by 
electronic glitches induced in the DISCSC data by gaps in other data types. The discovery in 
each case resulted from investigating an apparently significant (but spurious) negative lag. In ten 
cases the CCF did not consistently fit a peak in at least 50 of 101 iterations (the program found 
the "peak" at the edge of the CCF), and therefore the burst was eliminated. The final sample 
with measured lags contained 1437 long bursts. 

Figure 2a illustrates the computed lags versus peak flux in a linear-log plot. In this rendition, 
the lags (uncorrected for cosmology) for the BATSE bursts with redshifts (see Figure 1) span a 
relatively narrow zone - but still a factor of ~ 40 (0.01-0.37 s). Yet the longest lag for a burst 
source at high redshift is still 25 times shorter than the longest lag in the sample, ~ 10 s. But, 
relatively short lags dominate the BATSE sample (see Band 1997): more than 1000 bursts have 
Tiag < 350 ms. The Xiag-Fp plane is shown divided by solid lines into three regions containing: 
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120 bright bursts with Fp > 8 photons cm s and Xiag < 0.25 s, whose error bars are comparable 
to or smaller than the symbols in Figure 2a; 945 dimmer bursts with the maximum Xiag ranging 
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from 0.25 s to 0.45 s (as Fp decreases from 8 to 0.25 photons cm s ), where the error bars 
increase signiflcantly, with some extending to — 1 s; and 372 bursts with longer lags, and error 
bars typically Axiag/tiag ~ 25%. The point for GRB 980425 / SN 1998bw is circled (xiag = 2.8 s, 
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Fp= 0.9 photons cm s ). Figure 2b with the lag coordinate magnified shows more clearly the 
dense region of low peak flux and short lag. The dispersion in lag towards negative values is 
attributable completely to measurement error, as demonsfrated in section 2.3. 

2.2 Long-lag Bursts 

The solid line with negative slope in Figure 2b (dividing the second and third regions of the 
Xiag-Fp plane) was positioned to take into account the larger lag errors at lower peak fluxes. This 
is an attempt to separate the "long lag" bursts, which begin to dominate at low peak flux, from 

2 1 

the dim, short lag bursts. Notice that below Fp ~ 0.6 photons cm s the frequency of bursts of 
all lags clearly diminishes. This reflects the decreasing completeness of the BATSE sample as the 
trigger threshold is approached. More meaningful is the ratio of number of long-lag bursts (region 
3) to short-lag bursts (regions 1+2), shown in Figure 3 as a function of Fp in dyadic steps. This 
ratio increases dramatically from zero for the brightest bursts, to unity near trigger threshold. 
Figure 4 illusfrates the integral peak flux-frequency disfributions (Log[N>Fp]-Log[Fp]) for the 
long-lag bursts, short-lag bursts, and bright short-lag bursts. The long-lag bursts follow a -3/2 
power-law over ~ 1_ decades in Fp (to much lower in Fp than do the bright bursts), with the 
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inevitable rollover consistent with trigger threshold effects. The implications are that a relatively 
nearby GRB component of low luminosity is being detected, and that this component begins to 
dominate the frequency of burst detection near BATSE threshold. These long-lag bursts have the 
softest spectra of all three regions delineated in Figure 2a (Bonnell & Norris 2002). 

GRB 980425 is the canonical long-lag, soft-spectrum, ultra-low luminosity burst, its source 
lying only ~ 38 Mpc distant (Galama et al. 1998). Therefore, it is necessary to examine the sky 
distribution of the long-lag bursts in Supergalactic (SG) coordinates (Hudson 1993). Table 1 lists 
the quadrupole moment (Q) and error (8q) for long-lag bursts for xiag > 0.5-3.0 s in 0.5-s steps, 
along with number of bursts (Ngrbs) ^nd formal significance expressed in standard deviations (see 
Hartmann et al. 1996). The quadrupole moment for the 1065 bursts with short lags (those bursts 
in regions 1+2 discussed above) is Qshort lag = -0.005 □ 0.009. Since these bursts are believed to 
be at cosmological distances, Qshort lag constitutes an empirical measure of the quadrupole 
moment for BATSE exposure. In Table 1 the Q values were reduced by Qshort lag? ^nd the 
exposure error was propagated in quadrature to obtain the 8q values. For very long-lag bursts, 
Tiag > 1.5 s, the significance of the quadrupole moment is as high as 2-2.7 o before dropping as 
the number of bursts becomes small. Figure 5 shows the SG distribution of the 72 bursts with 
Tiag > 2 s; the concentration of this sample towards the SG plane is evident, with _ of the sources 
occupying that half of the sky between -30nand +30nin supergalactic latitude. The center of the 
Virgo Cluster is indicated by the large open circle. GRB 980425 / SN 1998bw hes at (bsc, Isg) = 
(20.8, -100.9). 

Taken together, the guilt by association with GRB 980425, the approximately -3/2 power 
law for all long-lag bursts, the preferential detection of these bursts near BATSE threshold, and 
the tendency of very long-lag bursts to follow the main feature of the nearby matter distribution 
- while not incontrovertible evidence - suggest that very long lag implies ultra-low luminosity. 
The implications of this conclusion are explored in the Discussion. Also indicated in Figure 5 is 
the position of GRB 971208, probably the BATSE burst with the longest spectral lag (visual 
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estimate ~ 20-30 s) and median peak flux, ~ 1.3 photons cm s . Its mono-pulse temporal 
profile spanned an Earth occupation (Connaughton et al. 1997), and was so long that it was not 
included in the lag analysis performed here since a usable background could not be fitted. 

2.3 Modeling Spectral Lag 
The next step is to model the Xiag-Fp scatter plot of Figure 2a with few parameters. The 
motivation is to realize a noise-free representation of the form N(Tiag, Fp) = f(Fp), which can be 
used to estimate redshifts and luminosities. Construction of a satisfactory representation of 
N(Tiag, Fp) is facilitated by dividing the Xiag-Fp plane into peak flux ranges in dyadic steps, as 
illusfrated in Figure 6. The increase in size of measurement error is evident from the larger 
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dispersion in lag near xiag = for successively dimmer Fp ranges. However, the dispersion 
cannot be completely reproduced by assuming that lags for dimmer bursts follow the same 
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distribution as lags for bright bursts (Fp > 8. photons cm s ) but with lower S/N. This is 
demonstrated in Figure 7, where the lags for the 120 brightest bursts have been recomputed for 
the three lowest Fp ranges (0.25-0.5, 0.5-1.0, and 1.0-2.0 photons cm"^ s"') but with peak 
intensities and S/N levels chosen randomly from bursts in the respective Fp ranges (see Norris et 
al. 1994 for a description of the Poisson S/N-equalization procedure). The bright-burst lag 
histograms are gray-filled and their peaks are normalized to the peaks of the solid-line histograms 
of the dimmer subsamples. In the core, Xiag < 0.5 s, the distributions for S/N-equalized bright 
bursts are significantly narrower than the distributions for the three dim Fp ranges, indicating that 
dimmer bursts tend to have a longer lag cutoff in the core component. Presumably, this reflects 
the detection of lower luminosity bursts at lower Fp. 

A model for the observed lags must reproduce the increasing fi-equency of longer lags at lower 
peak fluxes, both in the core and extending across the tail of the distribution to Xiag ~ 10 s. It is 
possible to reproduce both aspects using one power-law distribution, but with a long-lag cutoff, 
Tjiiax^ which increases with decreasing Fp. 

The minimum lag for bright bursts, XminO^ is another flt parameter. This minimum then would 
correspond to some maximum luminosity. However, just like other spectral-temporal measures - 
e.g., average peak-aligned proflle, duration, pulse interval - for which observed time-dilation 
trends with peak flux have been observed (Norris et al. 1994; Bonnell et al. 1997; Deng & 
Schaefer 1998), spectral lag most probably exhibits a similar average stretching trend as peak flux 
decreases: lag is a surrogate measure for pulse spectral evolution. Wider pulses have longer lags, 
and pulse width is affected by the twin cosmological effects of time dilation and spectral redshift. 
At this point in our understanding, it is not clear whether the observed time dilation is extrinsic 
(due to expansion of the Universe) or intrinsic (due to burst djoiamics and/or observer viewing 
angle). Yet the question is presently irrelevant, since the modeling should take into account the 
observed trend regardless of cause. The problem is how to represent an observed time dilation 
given that most spectral lags are short compared to most other timescales in bursts, and that lag 
measurements at low peak flux have large associated errors. The several measures of extrinsic 
time dilation/spectral redshift depend (nonlinearly) in different ways on the two effects (Norris 
1996), but yield comparable trends with peak flux. One reliable expedient is to use time-dilation 
factors (TDFs) measured for T50 durations; the uncertainties of this method are well understood 
(e.g., Bonnell et al.). Figure 8 shows the average T50 durations for 10 peak flux ranges, each 
group containing ~ 100 bursts. Only the bursts in regions 1 and 2 of Figure 2 were used for the 
purpose of estimating time dilation; the long-lag bursts were excluded based on the conclusions 
reached in section 2.2. A quadratic form was fitted to the TDF-Fp trend, and Xmin was modeled 
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as a (fixed) fimction of peak flux to yield an empirical stretch factor, Sobs = a + b Fp + c Fp , 
with a = 1.47, b = -0.722, and c = 0.277. 

Five parameters are then required to describe the modeled lag distribution - XminO^ "CmaxO^ P, 
and Fpo - such that 

N(Tniodel5 Fp) ~ "I-model 
with Tmin(Fp) = TminO Sobs Fp < 25 

"t-minCFp) ~ t-minO Fp [] 25 

and Tinax(Fp) = (Fp/Fpo)'' . (2) 

Integration and inversion results in the Monte Carlo formulation for the modeled lag, 

'''model ~ '''min { 1 ~R [1 ~ (t-max/^min) ] } ^ ^ (3) 

where R is a uniformly distributed random deviate. Thus, while a single power-law index, a, 
describes the lag distribution, the relative density of short lags thins out as XmaxCFp) increases. 

Several attempts were made to model the lag errors analj^ically as a function of lag and peak 
flux; none provided adequate results. Instead, the model errors were gotten fi-om the real error 
sample: For each model lag, an error for an associated real lag was chosen randomly fi-om a region 
in the Xiag-Fp plane within a factor of 01.20 in both coordinates. For a few cases (~ 10) the 
region was not populated and so the coordinate ranges were enlarged 10% per iteration until at 
least one point was included in the region. The "best" representation was then computed by 
minimizing a statistic similar to absolute value norms discussed by Scargle (1981), 

<d> = {2 IXmodel - "Creall / (I'^modell + l'Creall)/2}/Niags (4) 

searching for the minimiim over the 5-parameter model space. A minimization is not 
appropriate in this situation, since the lag errors are not Gaussianly distributed over the Xiag-Fp 
plane, but rather are strong fimctions of both peak flux and lag (and implicitly, a fimction of the 
disparate time profiles as well). The <d> statistic instead minimizes fi-actional distance in 
model-real lag pairs (the Fp coordinates being equal), where lag varies over ~ 3 orders of 
magnitude. The problem remaining in order to employ this statistic was, how to one-to-one 
associate a real lag and a model one? This problem was solved by decimation, finding the 
smallest distance in the Tiag-logio(Fp) plane for a pair, with the decimation starting at the furthest 
distance from the approximate centroid, {xiag, Fp} = {0.84, 1.2}. The decimation thus proceeded 
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from the sparsely populated outer regions to the cenfroid, elirninating the closest model-real lag 
pair per step. 

The resulting "best fit" set of model lags to the real lags in the Xiag-Fp plane was 



Tmodel (S) = 0.075 { 1 - R [1 - (WTminr-''] } 

with T„,ax (s) = 0. 125 (Fp/25.) 



-1/1.15 



(5) 



A typical realization of model lags is illusfrated in Figure 9a, as gray-fill histograms on top of the 
measured lags (solid histograms) for the same six peak fiux ranges shown in Figure 6. The general 
agreement appears satisfactory. Notably, the region near zero lag for 0.25 < Fp < 0.50 is 
adequately modeled. Without inclusion of the Sobs time-dilation factor, lags near zero tend to be 
overproduced in this lowest Fp range. Note that as Xmax increases with lower Fp - generating 
longer lags in the core - since fractional errors are large for xi^g ~ s, the negative lags at low Fp 
are fairly well reproduced by the model. In Figure 9b the same picture is illusfrated except 
without errors for the model lags. Now the effect of the fixed value for a is apparent in the 
persistent peak near positive, short lags. In Figure 10 the model lags with errors are plotted in 
the Tiag-Fp plane. Fidelity in representing the scatter plot for real lags (cf Figure 2a) appears 
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acceptable except for a handfiil of (real) long-lag outliers with Fp > 3 photons cm s , which are 
not reproduced due to the fixed dependence of Xmax on Fp - a minor defect for most purposes. 



The last step is to correct the modeled lags (without noise) for time dilation and specfral 
redshift. The cosmological fransformation can be realized by a procedure that generates a 
luminosity, Lo, from a lag-luminosity relation, where lag depends on z. The redshift is varied 
until the modeled peak fiux agrees with that observed (our Fp, V1-V2 : 50-300 keV). The tutorial 
on distance measures in cosmology by Hogg (2000) is instructive, in which Eqs. (14—16, and 22) 
relate differential peak flux to luminosity distance, redshift, and differential luminosity in the 
source frame: 



An estimate of the luminosity in the "deredshifted" bandpass, L(i + z) ^1^2, can be gotten from the 
ratio of fiux in the (1 + z)x(vi-V2) bandpass to a "bolomefric" fiux, multiplied by Lq. To 
calculate the flux ratio a simple model for the average GRB spectrum was assumed, a broken 
power law with fixed low and high energy specfral indices. The parameter values adopted are the 
approximate modes of the analysis by Preece et al. (2000) for the Band model (Band et al. 1993): 



2.4 Estimating Cosmological Corrections 



Fp,model ~ Fp^ ^1^2 - (1 + z) L(i + z) ^1^2 



I An Dl^ . 



(6) 
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N/E) = E"" = E"^ E<Ejoin = 230keV 

= E-P = E-2-25 EQEjoin. (7) 

The difference between a Band model and broken power-law is of no consequence in this 
treatment. In fact, the differences in results obtained for p = 2 or 2.25 are also negUgible. Then 
the expression for the band-limited luminosity in the source frame is 

(1 + z)v2 Emax 

L(i+z)vi^2 = Lo EN/E)dE / EN/E)dE (8) 

(1 + z)vl Emin 

where the limits on the y-ray bolometric flux were chosen to be Emax = 25 keV and Emin = 20 
MeV. While only ~ 20 burst spectra have been measured with the EGRET calorimeter, the 
average GRB spectrum is inferred to continue as a power law into the 10s of MeV regime 
(Dingus 2000). In principle, a full treatment should take into account spectral variation 
correlated with lag and peak flux. But since the variation in high energy power-law index noted 
above influenced the results negligibly, this approach with a fixed spectrum is probably a 
reasonable approximation. 

The lag-luminosity relation needs to be augmented in this model to include lags longer than a 
few 100 ms, i.e., to include the long-lag bursts of apparently ultra-low luminosity. This other 
branch of the "GRB HR diagram" (NMB) is empirically constrained only by GRB 980425 and 
the absence of any bursts on the main branch with lags longer than ~ 300 ms. Hence, a partially 
arbitrary conjecture to satisfy this lack of constraint would have the hj^othetical low-luminosity 
branch commence near 350 ms, or L53 « 2.2 x 10 (see also the similar, but theoretically based 
conjecture of Sahnonson 2001), and thus 

Lo = L53= 1.3 X [Tiag/0.01 s]"^-^^ 0.003 < Tiag < 0.35 

= 7.8 X [Tiag/0.1 s] ""^-^ 0.35 < Tiag . (9) 

The second defining point for the low-luminosity branch is GRB 980425, for which L53 ~ 1.25 x 
10"*^ ergs s"^ and the measured Xiag = 2.8 s (cf. NMB, where the lag for the whole profile was 
estimated visually as ~ 4.5 s). 

The remaining ingredient is correction of the model lags for the extrinsic effects of time 
dilation and spectral redshift. The former correction factor is trivial, (1 + z). The correction for 
redshift of temporal structure from the fiducial 25-50 keV and 100-300 keV bands in the source 
frame, to lower energies in the detector frame, can be only approximately modeled using 
measurements of the brightest BATSE bursts' average temporal structure (Fenimore et al. 1995; 
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Norris et al. 1996). These bursts will tend to be at lower redshift than dimmer, short-lag bursts. 
At a redshift of 1 + z = 5 (10), the 25-50 keV band is shifted to 5-30 keV (2.5-15 keV). No 
quantitative narrow-band measurements of the relative widths of temporal structure have been 
reported in the X-ray bands. Visual inspection by this author of BeppoSAX and older X-ray 
temporal profiles of GRBs appears to indicate that pulse structures broaden faster per 
logarithmic interval at these low energies than at BATSE energies. Here for expediency, the 
average pulse width of bright BATSE bursts is utilized to correct roughly for redshift, W<puise> D 
g-033 (-[sjQj.j.js Qi 1996), but with the exponent modified by a sigmoid function at high redshifts. 
Then the combined cosmological correction applied to Xmodei of Eq. (5) is 

Smodel=(l+z)-\l+z)'^ (10) 

with K = 1 - (2/3){2/[l + exp(z/zijiii)}. The expression reduces to k = 1/3 at z = 0, and 
asymptotes to unity as z surpasses ziim, where it nearly equals and counteracts time dilation (the 
effect of spectral redshift could be even larger: k could grow larger than 1). The lack of 
knowledge of k at low energies represents the largest source of uncertainty in the whole 
treatment. Hence the model is evaluated for two illustrative trial values of ziim, 5 and 20. 

The iteration procedure then proceeds as follows: Assume the cosmology {Qm, ^ {^-3, 
0.7} with Ho = 72 km/s/Mpc (Freedman et al. 2001). Guess an initial redshift (zq), compute the 
cosmology-corrected lag fi-om Eq. (10), and integrate the luminosity distance. The source fi^me 
bandpass and lag-predicted luminosities defined by Eqs. (8) and (9) yield a model peak flux (Eq. 
[6]) in the detector fi-ame, Fp modeh to be compared with Fp obs- The next guess is generated fi-om 
Znext = Zoid (Fp,modei/Fp,obs)~- The procedure converges (1% accuracy in Fp) within a few 
iterations yielding the same Zfmai whether Zq = 1 or 0. 

Figures 1 1 through 13 illustrate the immediate results obtained for spectral lags, redshifl:s, and 
proper distances, respectively, for ziim = 5 (solid histograms), and 20 (dashed histograms). Note 
that these modeled distributions are representations for the BATSE sensitivity-limited sample of 
long bursts. A common feature in these figures is the necessary convergence of the distributions 
for both values of ziim for the putative nearby sources of long-lag GRBs. Also, all the 
distributions peak near values corresponding to large redshifts. The cosmology-corrected lags 
(Figure 11) peak near 7 ms (zHm = 20) and 13 ms (ziim = 5), the majority occurring in either case 
at Tiag < 50 ms, implying that the observed luminosity distribution is peaked in the neighborhood 
of L53 ~ 1. Figure 12 shows the redshift distribution plotted as (1 + z), thereby artificially 
emphasizing the nearby sources. For ziim = 5 (20) the sources at cosmological distances are most 
numerous near z ~ 10 (20). This uncertainty in where the GRB rate-density peaks just affirms 
our lack of knowledge of the appropriate form of the spectral correction for high redshift. Figure 
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13 again emphasizes the preponderance of sources at large proper distance. But the other 
interesting feature is that ~ 90 sources are modeled to He at distances < 100 Mpc. A more 
accurate estimate of their distances would require better definition of the steeper lag-luminosity 
branch defined in Eq. (9). These apparently nearby sources may have interesting ramifications 
for other astrophysics topics such as gravitational radiation fi-om aspherically collapsing objects 
and ultra-high energy cosmic rays, and for the future Swift mission (see Discussion). 

In Figure 14 the scatter plots of peak flux versus redshift for the two values of have the 
expected appearance: The hard limit of Lgutoff is manifest in the upper right boundary of each 
populated region. For = 20 the sources actually pile up at this boundary, reflecting the 
increase of redshift to accommodate short-lag, high-luminosity bursts, suggesting that with ziim = 
20 the (1 + zf factor in Eq. (10) ramps up too slowly. For ziim = 5 (20) the median z increases 
from ~ 2 (3) for bright bursts to ~ 10 (20) for the dimmest bursts. While star formation rates are 
highly uncertain at these redshifts, investigators do not currently envision that star formation 
peaks at z > 20 (e.g., Madau, Delia Valle, & Panagia 1998), and so a ziim as high as 20 in Eq. (5) 
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is not favored. Below Fp ~ 4 photons cm s , the concentration of low-luminosity, long-lag 
bursts appears near z = 0. Besides these obviously expected features, there is a sparsely 
populated triangular region at low peak flux between the nearby low-luminosity GRBs and the 
high-luminosity cosmological GRBs. This sparsely populated zone (z < 2.3 near BATSE 
threshold) is consistent with a "seeing through" of the GRB population since lower redshift 
sources should continue to be detected at the same peak flux where higher redshift sources are 
represented. 

It is interesting to compare the current distribution of 24 redshifts obtained from GRB 
afterglows or host galaxies, shown in Figure 15 (see the GCN alerts, maintained by Scott 
Barthelmy: http://gcn.gsfc.nasa.gov/gcn/gcn3_archive.html), with the Fp-z scatter plot in Figure 
14b. The mode and median for determined redshifts continue to hold near z « 1, with ~ 67% of 
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the sample contained within 0.4 < z < 2; whereas for the model above Fp ~ 10 photons cm s 
the median is z ~ 2.5. Approximately 120 bursts have been detected by BeppoSAX, RXTE, 
HETE-II, and/or the interplanetary network since the "afterglow era" began in February 1997. In 
the overwhelming majority of cases when X-ray observations were made. X-ray afterglows were 
detected, but for only ~ 40-50% of these were optical afterglows detected, perhaps half the time 
leading to host galaxy identifications with redshifts associated to the GRB (Kevin Hurley, private 
communications). Thus obscuration at the source or redshift of the spectrum have prevented 
optical afterglow detection, probably preferentially for higher redshift GRB sources. 
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3. COMPARISON OF MODEL AND THEORETICAL LUMINOSITY DISTRIBUTIONS 



For the BATSE sensitivity-limited sample of long duration (T90 > 2 s) bursts, the luminosity 
distribution (znm = 5) derived from the model lags is illustrated in Figure 16. On the right side of 
the figure the seven cascading dotted/solid histograms are for the redshift-limited cuts: z < 30, 
20, 10, 5, 3, 2, and 1. The plot is rendered in equal log-spaced intervals; hence the slopes of the 
two fitted power laws that are illustrated appear one unit flatter than for linear-spaced intervals. 
The distributions for z < 5, 10, 20, and 30 show increasing numbers of bursts at higher 
luminosities - the observed luminosity distribution peaks near L53 ~ 1. But this just reflects the 
deflcit of detectable lower luminosity bursts at higher redshift. GRBs with L53 < 10~^'^ 
(leftmost dotted histogram, z < 0.024, d < 100 Mpc) comprise the subsample with xiag < 2 s. 
Recall that the model employs the two-branch lag-luminosity relation expressed in Eq. (9), which 
breaks at Xiag = 0.35 s, or L53 = 2.2 x 10 . The number of long-lag bursts per logarithmic 
interval decreases as lag increases across the break point. The trend can be roughly described as a 
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power-law, dN/dxiag ~Xiag ' over 0.01 < Xiag (s) < 3 (see Figure 11). Thus it is mostly the 
steeper power law of the low-luminosity branch which gives rise to lower fi-equency of GRBs 
per logarithmic interval for L53 < 2 x 10 . 

As noted, the model luminosity distribution reflects the sensitivity-limited BATSE sample. 
However, since we know something about GRB redshifts, a high-luminosity volume-limited 
regime can be rendered fi-om this picture for particular cuts in luminosity and redshifl:. From 
inspection of Figure 14 the suggestion (section 2.3) was that BATSE has sampled through the 
GRB population for z < 2.3. Consider a source number density ri(L)dL = dN/dV, where ti(L) 
may implicitly include a luminosity dependence for viewing angle, a profiled jet, etc. The number 
of sources actually detected (in Euclidean space) is then 

p 

N(L)dL = [}i(L)dL 4jr ? dr (1 1) 



3 3/2 

where p = min{Ro , [L MjtFp thres] } and Fp thres is the peak flux at BATSE trigger threshold. 

1/2 

Thus for Rq < [L MjtFp thres] ■> a sample is volume limited, N(L)dL r|(L)dL, and we see the 
intrinsic luminosity distribution. (Any evolution with redshifl: is ignored in evaluating the 
average N(L); evolution may be a discriminant for GRB models, but it is not examined here). 
Beyond z ~ 0.2, the luminosity distance must be used since it significantly exceeds the Euclidean 
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distance. From Eq. (7) the average photon energy is 135 keV, and so Fp thres " 5.4 x 10 
ergcm"^ s"\ Thenforz<2 andL53 > 1.5 x lO"'^, Dl(z=2) =15.1 Gpc, and [L /43tFp thres] = 
15.2 Gpc, and so this constitutes a volume-limited regime, integration of Eq. (1 1) yields N(L)dL 
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n ri(L)dL. Over the volume-limited regime the distribution has a fitted power-law slope of -0.82 
(solid straight line), so 

dNvoi/dL[] L"^•^ and 

Ti(L)nL"^-^ 10"^-^<L53< 10°-^ (12) 

—3 1/2 

Whereas, for z < 1 and L53 < 2.6 x 10 , Dl(z =1) = 6.4 Gpc and [L MjtFp thres] = 6.3, and 

this regime is sensitivity limited - the sources are progressively undersampled as luminosity 

3/2 

decreases across this regime. Integration of Eq. (11) now yields N(L)dL r|(L)L dL. For the 
sensitivity-limited regime the fitted power-law slope is +0.015 (dashed straight line), so 

dNsen/dLDL"^ and 
Ti(L) D L"^^^ 10"^-^ < L53 < 10"^-^ . (13) 

(If the lag breakpoint for the two-branch lag-luminosity relation in Eq. [9] is increased from 0.35 

s to 0.6 s, then L [xjag/O.l s]~^'^, and the fitted power law slope becomes -0.015; hence the 

slope for the sensitivity-limited regime is not very dependent on the precise break point 

position.) Below L53 ~ 10"^'"* the detected source frequency falls off, indicating either that such 

low-luminosity GRBs are below the BATSE trigger threshold, or perhaps that the intrinsic 

distribution cuts off. For reference, the burst with lowest determined redshift and luminosity, ~ 

1.25 X 10"^L53 for GRB 980425 (z = 0.0085, d = 38 Mpc, Fp = 0.9 photons cm"^ s'\ would 
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have been detected by BATSE at the sample threshold, Fp ~ 0.25 photons cm s , at ~ 72 Mpc. 

2 1 

Recall the handful (~ 15) of real long-lag outliers (Fp obs > 3 photons cm s ) not reproduced 
at comparable Fp model by the model lags (cf Figure 2a and Figure 10) due to the fixed (versus 
fuzzy) dependence of Xmax on Fp. This represents a truly minor defect in the modeling among 
the ~ 210 bursts with L53 < 10 over 5 decades in luminosity: According to their real positions 
in the Xiag-Fp plane, about half would tend to fall in the underpopulated region < 4 x 10 L53, 

_3 

and the rest would be disfributed fairly uniformly up to ~ 10 L53. 

3 . 1 Luminosity Distribution via Range of Viewing Angle, Constant F 
Eqs. (12) and (13) can be compared to expectations from purely relativistic kinematic 
proposals that viewing angle. By, off the jet axis is the governing factor for the observed 
luminosity distribution (see Sahnonson 2000; 2001). For instance, in Salmonson & Galama 
(2001) the luminosity scales (neglecting redshift) as the perceived Doppler factor, L D ~ 2r / 
(1 + [(By - Bjet)r] ) for By « 1 and F » 1, i.e., near the jet cone, Bjet- But for By large 
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compared to Bjet (~ few degrees), L ~ D . and the GRB appears as a low-luminosity source. In 
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the latter case, for F ~ few and constant, L ~ [SyF] 6v . The perceived source density 
distribution in the pure viewing angle scenario is Ti(L)dL [] sin(6v)d6v. Even for 6v < 40Q 
ri(L)dL ~ eydOv, and so Ti(L)dL ~ L'^^'^dL, and N(L)dL Q Ti(L)L^^^dL ~ L+^^^dL. This slope, 
more than one unit steeper than Eq. (13), would result in a factor of ~ 3000 overproduction of 
GRBs at 10 L53 compared to 10 L53. This is consistent with the approximation made in the 
solid angle integration, given that By varying over a factor of 10 (say, 6v,min ~ 40 > 6jet) would 
result in a dynamic range of ~ 10^ in luminosity. Thus the extended flat distribution provides a 
good discriminant in the sensitivity-limited regime against the purely relativistic kinematic 
explanation for the low end of the modeled GRB luminosity distribution. 

For the high-luminosity, volume-limited regime, 10"^'^ < L53 < 10^'^ with z < 2, the picture is 
less straightforward to interpret, since a significant fi-action of bursts will still be observed within 
the jet cone - and if the beaming fi-action varies appreciably (possibly correlated with F), we are 
in mixed territory (which may be the actual situation). Still, keeping with the purely relativistic 
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kinematic interpretation and 6jet< few degrees, since F» 1, L~ [(6y - 6jet)F] = [6aF] . Then 
via the same reasoning as above, ri(L)dL ~ L dL, not significantly different fi-om Eq. (12). 
However, any realistic dependence of the Lorentz factor on angle with respect to the jet axis (F 
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6jet , X > 0) will render r|(L) steeper than L . So again, the purely relativistic kinematic 
explanation is not favored to explain the distribution of highly luminous GRBs. 

3.2 Luminosity Distribution via Variable Beaming Fraction 
In the pure beaming fi-action scenario, the luminosity and energy radiated into a differential 
solid angle are taken to be constant across the face of the jet cone, and negHgible outside it. Also, 
the luminosity radiated into 4n steradians is found to be an approximate invariant (Frail et al. 
2001), such that D fb D S2jet/43t. To reproduce the modeled r|(L) dependences in Eqs. (12) 
and (13) requires distributions in opening angle, or jet-cone solid angle, dN(Qjet)/dQjet D ^^jet^- 
It follows that A, = 0.2 for the volume-limited regime, 

dN(Qjet)/dQjet D S2jet"" ' , 6 < 200 (14a) 

and X = -0.5 for the sensitivity-limited regime, 

dN(Qjet)/dQjet D ^2jet^"•^ e>20n. (14b) 

The break at 6 = 20n(fb = 0.06) is evaluated for xiag = 0.35 in Eq. (1), corresponding to the break 
in the assumed two-branch lag-luminosity relation of Eq. (9). For the high-luminosity GRBs, Eq. 
(14a) predicts a slowly decreasing distribution for Qjet. Thus, as the opening angle increases 
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from the rairiimvim to maximum values derived by Frail et al. (2001) and PK (6jet ~ 2-20Q, 
dN(Qjet)/dQjet decreases by only a factor of ~ 0.6. For the low-luminosity GRBs, Eq. (14b) 
predicts that the number of bursts actually increases approximately linearly with 6jet (since 

Qjet/43T « _ Ojet^). 

Now, employing the two branches of Eq. (1) that relate spectral lag and beaming fraction (and 
remembering that the relation for the long-lag branch is merely a toy model), the model lags of 
Figure 1 1 can be translated into numbers of bursts which would have been detected by BATSE if 
we had observed them all within their jet cones. This distribution, N[fb(Tiag)] / fb, is plotted in 

Figure 17 for ziim = 5. The implied number of highly luminous bursts (fb _ 3 x 10 , Xiag < 0.067 
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s, L > 10 ' L53) - most of them sufficiently luminous to be seen across the Universe - would 
need to be reproduced by a successfiil physical model. Only ~ 72% (1437/~2000 = f^se) of 
BATSE long bursts were analyzed in this study. The corrected rate is Ngrb _ 6.6 x 10^ / 
(9.2xfBATSEfuse) 5^ ^ where the average live fraction for BATSE was fsAxsE ~ 0-48 (Hakkila et al. 
2002). Therefore the rate of highly luminous bursts, irrespective of actual jet axis orientation, is 
estimated to be Ngrb > 2 x 10^ 50-"^ 

3.3 Luminosity Distribution via a Profiled Jet 
Models invoking a profiled jet with F, and therefore luminosity, decreasing off the jet axis 
have been developed recently by several authors. This dynamical mechanism can account for the 
anti-correlation between afterglow break times and luminosity, and possibly spectral evolution in 
the y-ray phase (MacFadyen, Woosley, & Heger 2001; Sahnonson & Galama 2001; Rossi, 
Lazziti, & Rees 2002; Zhang & Meszaros 2002). Consider Ljet D with y > 0. Then the 
solid angle dependence r|(L)dL ~ 0yd0 y yields Ti(L)dL ~ @y^^6L, and Ti(L)dL ~ L~^~^^^dL. 
Agreement with the model distribution in the high-luminosity regime of Eq. (12) requires -1-2/y 
= -1.8,orY = 5/2: 

LjetDev"'^'. (15) 

Within the profiled jet cone the observed luminosity is integrated over the observer viewable 
angle ~ and so the luminosity dependence is an integration of dependences on r(0y) as well 
as on 0v, which are implicit in Eq. (15). 

To estimate a lower limit on the number of highly luminous GRBs, a working model for the 
profiled jet must be assumed. For specificity, and to include the highest luminosity burst 

—5/2 

produced by the model lags, take L(0v) = 3.6 L53 x [0y/0o] and two values for the nainimvim 
jet cone radius, 0o = 1.50 and 3[\ (0.026 and 0.052 radians, respectively), below which the 
luminosity remains constant (to avoid divergence) at L = L(0o). Jets not oriented toward the 
observer increase the number of GRBs preceived to have L(0y) by a factor 43T/[2[]2jr0vd0y], 
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since the fraction of bursts with L > L^^v is 43t/[2Q^^^'^MQv] (the extra factor of 2 accounts for jets 
being two-sided). Figure 18a shows the resulting distribution of GRBs within the BATSE 
sensitivity reach for the cutoff L > 10"^^L53, irrespective of jet axis orientation with respect the 
observer. If the source density were slightly steeper than expressed in Eq. (12), say r|(L) L , 
then Y = 2. Figure 18b shows the distributions for Ljet D ®v and the same two values of Oq. 
GRBs with L > 10"'^L53 are volume-limited within z < 2, while the most highly luminous GRBs, 
up to L5 3 ~ 1, are observable even to z ~ 20. To compare with the rate calculated above for 
beaming fraction with L > 10 ' L53, for Ljet D ®v the total numbers are 2.7 x 10 (7 x 10 ) 
for 00 = 1.5n(3Q. For Ljet D ©v the corresponding totals are 2.1 x 10 (5.4 x 10 ). Then the 
rates are Rgrb > Ntot / (9.2xfBATSEfuse) F"^ « 0.66-0.85 x 10^ yr"^ (0o=1.5Q and 1.7-2.2 x 10^ 
yr~^ (00=3 Q. The lower limits on rates for the profiled jet and beaming fraction scenarios are 
compared with estimates of the SN Ib/c rate in the Discussion. 

For the sensitivity-limited regime, Eq. (13) implies Ljet D , a flatter dependence than 

for the highly luminous GRBs. Assuming that T ~ few for the long-lag GRBs, the observable 
y-ray cone (9^ ~ I^^) can be substantial fraction of lit. Hence for the low-luminosity GRBs, a 
large fraction of jets with axes not directly oriented toward the observer are still observed. 

4. DISCUSSION 

In the pre-afterglow GRB era it was often remarked that burst temporal-specfral behavior 
was chaotic - unpredictable - zoo-like, and discussions were held at workshops to attempt 
useful classification schemes. But prior to BATSE, such classifications were uninterpretatable 
because GRBs were largely unpredictable - on their duration timescales. However, the 
sensitivity of the BATSE experiment has allowed us to see and map an important organizing 
principle: pulse spectral evolution. Within a given burst the individual pulses tend to have 
comparable widths and spectral lag. But in the whole burst population, pulse widths and lags 
vary over a large dynamic range, > 10 . This variation appears to be inversely related to 
luminosity in long (T90 > 2 s) BATSE bursts. (A second organizing principle in y-ray time 
profiles, reported by Ramirez-Ruiz & Merloni [2001], is yet to be completely exploited: Time 
intervals near background level, between major emission episodes, are well correlated with the 
time interval of the following episode, sfrongly suggesting that a metastable accretion mechanism 
is fueling the burst process. Almost certainly, this correlation eliminates the external shock 
hypothesis as the energy release mechanism for the y-ray phase [E. Ramirez-Ruiz, private 
communication] .) 

The wide dynamic range in perceived GRB luminosity has been revealed by afterglow 
observations which yield associated redshifts. And now the afterglows (or their absence) also 
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manifest variable behavior - varying decay timescales and power-law break times. Many 
investigators have proposed unifying models which can explain the break times, dynamic ranges 
in luminosity and y-ray energy, and in some cases the y-ray pulse behavior (Frail et al. 2001; 
Panaitescu & Kumar 2001; Wang & Wheeler 1998; Hoflich, Wheeler, & Wang 1999; Nakamura 
1999; Woosley, Eastman, & Schmidt 1999; Salmonson 2001; loka & Nakamura 2001; Sahnonson 
& Galama 2001; Rossi, Lazziti, & Rees 2002; Zhang & Meszaros 2002). The economy of each 
of the three model classes - variable beaming fraction, pure relativistic kinematics (viewing angle), 
and profiled jets with graded Lorentz factor - is difficult to evaluate: each affords a reasonable 
explanation of the dynamic ranges in GRB energy and luminosity, such that the actual energy 
release into 4jr steradians is approximately constant, at least for GRBs at cosmological distances. 
In fact, reality may be a combination of the proposed effects. However, assuming a lag- 
luminosity relation, one can derive distributions for the large BATSE long-burst sample which 
constrain the three pure mechanisms described above. Here, I have analyzed spectral lags for 
1437 long (T90 < 2 s) BATSE bursts to near the trigger threshold, and modeled the lags to yield a 
noise-free representation. The model lags with cosmological corrections were used to compute 
distributions for redshift, proper distance, and luminosity. The main results, which may have 
further implications for GRB populations, are summarized: 

n A GRB subsample is identified - those with few, wide pulses, lags of a few tenths to several 
seconds, and soft spectra. The proportion of these long-lag bursts (Figure 3) increases from 
negHgible among bright BATSE bursts to ~ 50% at trigger threshold. Their integral 
size-frequency distribution (Figure 4) follows a -3/2 power-law over 1_ decades in Fp. 
Bursts with very long lags, ~ 1-2 < xjag (s) < 10, show a tendency to concentrate near the 
Supergalactic Plane (Figure 5) with a quadrupole moment of — 0.10 □ 0.04 (Table 1). GRB 
980425 (SN 1998bw) is a member of this subsample of ~ 90 bursts with estimated distances 
< 100 Mpc (Figure 13) and ultra-low luminosities (Figure 16). 

nin two luminosity regimes, the model lags yield approximate power-law scaling relations, 
Eqs. (12) and (13), for their respective distributions. These empirical relations are compared 
with expectations for jet scenarios which invoke variable beaming fraction with 
approximately constant total energy, relativistic kinematics (pure viewing angle effect), or 
profiled jets with variable luminosity (i.e., variable Lorentz factor) off the jet axis. In the 
high-luminosity volume-limited regime, 10^'^ < L53 < 10*^'^ with z < 2, the model 
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distribution is dNyoi/dL ~ L ' . For the low-luminosity GRBs in the sensitivity-limited 
regime 10"^'^ < L53 < 10"^'^, the model distribution scales as dNsen/dL~ L~\ 
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Both dependences can be reproduced with profiled jets, where Ljet D ©v for the high- 
luminosity bursts, and Ljet D in the low-luminosity regime. The beaming fi-action 
scenario can also fit the high-luminosity bursts with a fairly flat distribution for the jet-cone 
solid angle, dN(Qjet)/dQjet D ^^jet^'^- For the low-luminosity bursts a distribution that 
actually increases approximately linearly with 9jet (S2jet/43T « _ 9jet ) is required, 
dN(Qjet)/dQjet D ^^jet The pure viewing angle scenario produces a distribution which is 
too steep for the low-luminosity GRBs, dNsen(Ly dL L^^^^. For the high-luminosity 
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GRBs, dNvoi/dL L , only slightly steeper than the model distribution (0 L ' ). 
However, this would leave no room for inclusion of a profiled Lorentz factor, expected from 
modeling coUapsars (MacFadyen, Woosley, & Heger 2001), which would steepen the 
distribution. 

n The lag model of Eq. (4) and iteration procedure of section 2.4 for estimating redshift could be 
used in other programs to obtain results, e.g. for a volume-limited treatment of the GRB 
universe. Currently the largest uncertainty in any such treatment is our limited knowledge of 
spectral redshift corrections for high z bursts. In fact, in Eq. (9) ziim could be lower than 5, 
and K could be larger than unity for large redshifts. Therefore, it is prudent to view the 
results for ziim = 5 as much more useful for prediction compared to ziim = 20 where the 
extrapolation range is large. However, the luminosity distribution for high-luminosity bursts 
(Eq. [12]) should be unaffected by these uncertainties since it was constructed for burst 
sources with z < 2 and L53 > 10"^'^, a volume-limited regime for which the estimated redshift 
correction is reasonable. The modeled redshift distribution for GRBs peaks near z ~ 10, a 
result which has a large associated error, due to the uncertainty in correction for the effect on 
time profiles of redshifted GRB spectra beyond z ~ few. 

4.1 GRBs, SNe, Gravitational Waves, and UHECRs 
The GRB-SN connection began with GRB 980425 / SN 1998bw. The large number of long- 
lag bursts - like GRB 980425 in several respects - that are found mostly at low BATSE peak 
flux, suggest some requirements on SN Ib/c and coUapsar models (MacFadyen & Woosley 1999). 
Many of these long-lag bursts are implied by the two-branch lag-luminosity relation to be so 
close (within 10s of Mpc) that we should expect to identify frequently the responsible 
coUapsar/hypemova - regardless of what SN type it turns out to be - especially when alerted 
with Swift localizations (Barthelmy 2001). If these are SN Ib/c, and their predicted low Lorentz 
factors (r ~ 2-5: Kulkami et al. 1998; Woosley & MacFadyen 1999; Sahnonson 2001) are at 
least partially responsible for the ultra-low luminosities, then approximately all SN Ib/c are 
needed to account for long-lag bursts: Bloom et al. (1998) estimate ~ 0.3 SN Ib/c day~^ within 
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100 hg 5 Mpc. Over the ~ 9.2 years of CGRO operations, the BATSE sample used here would 
have had exposure to ~ 10 x fBATSfifuse ~ 350 SN Ib/c within the same volume. There are 372 
bursts in region 3 of Figure 2 with lags > 0.35 s, and 155 (72) bursts have lags > 1 (2) s. The 
sources of approximately 90 of these bursts are predicted (Figure 13) to lie within 100 Mpc. For 
r ~ few, the cone of observable y-ray emission is large, and so in combination with the flatter 
predicted dependence for a profiled jet (Ljet D ®v~^'^) than for the high-luminosity GRBs, the 
observer would be within the y-ray emission cone of most low-luminosity GRBs. Thus, a large 
fi-action, at least 25%, of SN Ib/c would participate in making long-lag GRBs, or else some 
unobserved nearby source type would be required. Notice that in a whole-Universe, volume- 
limited sample, the ultra-low luminosity bursts would therefore constitute the major variety of 
GRB (see discussion in MacFadyen, Woosley, & Heger 2001). 

Confirmation that the sources of long-lag GRBs are relatively nearby could come fi-om several 
directions. First, untriggered BATSE bursts - most with peak fluxes lying below the sample 
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threshold of this study (0.25 photons cm s ) - should be predominantly long-lag bursts (see 
Figure 3). However, the localizations for these dim events have large error regions, and their 
idenflcation as bona flde GRBs becomes more probabilistic with decreasing peak flux. Another 
possibility is that, since the matter distribution within 100 Mpc is not planar, cross-correlation 
of the positions of nearby galaxies (Hudson 1993) and long-lag GRBs may yield a higher 
confldence result than the simple quadrupole test. Also, the light curves of unusual SNe for 
which more than one detection exist may be extrapolated to the SN onset; however, this 
approach may be least likely to yield conflrmation since most SNe within 100 Mpc during the 
CGRO era went undetected (and still do), and such statistical treatments usually require ad hoc 
assumptions. The highest guaranteed approach will be real-time alerts provided by Swift, 
allowing detection of SNe if they are indeed associated with long-lag GRBs. 

Now consider the possible association of GRBs at cosmological distances with SNe. From 
van den Bergh & Tammann (1991) and Cappellaro et al. (1999) one may roughly estimate that ~ 
10% of all local SNe are type lb (or Ic). Whereas at high redshifts perhaps 20% of SNe are type 
Ib/c since these core collapse events were more frequent in the early Universe (see Figure 2 of 
Madau et al. 1998). Assuming ~ 1 SN s~Mn the observable Universe for all types, then ~ 6 x 
10 SN Ib/c occurred (at our epoch) during the CGRO mission. The pure beaming fraction 
scenario (Figure 17, section 3.2) predicts a rate > 2 x 10^ yr~^ for highly luminous GRBs (L > 
10"°-^L53), or 1.8 X 10^ GRBs during CGRO. Thus at least 1/30 of the SN Ib/c population at 
high redshifts - presumably the highly unusual and energetic ones - would be required to produce 
the most luminous GRBs. However, if the volume-limited luminosity distribution continues to 
scale as expressed in Eq. (12) beyond z = 2, then the total rate of GRBs with 10~^'^< L53 < 
10^'^, is ~ 4—5 times larger, or ~ 10^ yr ^ Then ~ 1/7 of the SN Ib/c population would be 
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required to produce the cosmologically distant GRBs in the pure beaming fraction scenario. 
Similarly, for the profiled jet scenario, the total GRB rates would be ~ 3 x 10^ yr~^ and ~ 10^ 
yr~^ for ininimviin jet cone radii of 0o=1.5nand 0o=3D respectively (Figure 18, section 3.3). 
Again, a large fraction of the SN Ib/c population would have to participate in making highly 
luminous GRBs in the profiled jet scenario. 

Results similar to those found here for the GRB redshift disfribution have been inferred by 
Fenimore & Ramirez-Ruiz (1999). In addition to the GRB rate-density peaking at large redshift 
(z ~ 10), Lloyd-Ronning, Fryer, & Ramirez-Ruiz (2001) find that the GRB luminosity 
disfribution scales as (1 + z)^ They conclude that the star formation rate continues to 
increase to redshifts at least as high as z ~ 10, in order to provide a sufficiently high progenitor 
rate for GRBs. Inferences on star formation in the early Universe based on GRB analyses may 
be compared with quasar observations in which the Gunn-Peterson frough is detected - evidence 
of the epoch of reionization. Djorgovski et al. (2001) and Becker et al. (2001) conclude that the 
reionization era was complete near z ~ 5-6. Simulations suggest that the reionization era spanned 
z~ 15toz~5, roughly consistent with GRBs appearing at redshifts of at least z ~ 10 (Gnedin 
2000; Ciardi et al. 2000; Umemura, Nakamoto, & Susa 2001). 

With the optimization of LIGO II, S/N levels of order few x 10"^"* Hz~- should be detectable, 
depending upon particulars of instrument enhancement (Fritchel 2001). For some tuned 
configurations the LIGO II sensitivity peaks in the ~ 100-500 Hz range. Modeling of 
axisjmmietric core collapse by Zwerger & Miiller (1997) using two-dimensional Newtonian 
hydrodynamic code predicts peak sfrains of order 10"^^ in the 100-800 Hz band for distances of 
10 Mpc. These expectations are reinforced by general relativistic modeling performed by 
Dimmehneier, Font, & Miiller (2001) which agree with amplitudes from Newtonian models 
within ~ 30%. Also, Fryer, Holz, & Hughes (2001) expect sfrains from bar modes of few x 10"^^ 
to X 10"^^ at 10 Mpc. For sources at 100 Mpc, the sfrains would be 10"^'*-10"^^ For a signal 
persisting in this band for a few cycles, the LIGO II sensitivity franslates into a sfrain sensitivity 
of ~ 3 X lO"^"* (300/few) - ~ 3 x 10"^^ for a S/N ratio of ~ 1:1. Thus a small fraction of the ~ 
90/(9.2fusefBATSE) ~ 30 long-lag bursts yr~^ (over several years) within 100 Mpc expected from 
the BATSE sample analysis may be detectable by LIGO II, at least if the signals are co-added 
(see Finn 2001). Since the long-lag GRB sources appear to be nearby, their signals would not be 
time-dilated beyond LIGO's sensitivity range. Detection will require computation of accurate 
chirp templates for core collapse events, and knowledge of the GRB occurrence time and location 
(an argument for maintaining all-sky monitoring of GRBs). 

The long-lag bursts might be sources of ulfra-high energy cosmic rays (UHECRs). Nagano & 
Watson (2000) provide an complete summary of the field. Proposals for GRBs at cosmological 
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distances (F ~ 10 -10 ) as the source of UHECRs have been made by Vietri (1995) and Waxman 
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(1995). Vietri shows that the cosmic ray energy flux > 3 x 10 eV expected from cosmological 
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GRBs is comparable to that observed, 5x10 ergs cm s sr . However, Stecker (2000) and 
Scully & Stecker (2002) argue that energy losses suffered by GRBs at cosmological distances 
would produce too sharp a cutoff in the UHECR spectrum, and would fall short in producing the 
flux by a factor of ~ 100-1000. The latest results from AGASA indicate a pronounced flattening 
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of the spectrum > 10 eV compared to results from previous experiments (Watson 2001). Thus 
the salient question is, how does the UHECR flux from ~ 30 GRB yr~^ within 100 Mpc and F ~ 
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few compare with that for 900 GRB yr Universe beamed at us with T ~ 10 -10 , and could 
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the closer GRB sources provide the enhancement seen by AGASA > 10 eV? There is the 
perennial problem of (an)isotropy. The AGASA results possibly indicate point origins, 
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implying nearby sources, but the interpretation is still open. For energies > 10 eV, source 
distances < 100 Mpc, and magnetic field scale lengths X ~ 10 Mpc, the arrival directions are 
scrambled if B > 10 G, a typical exfragalactic field sfrength believed by many investigators (see 
Kronberg 1994, on extragalactic magnetic fields). Farrar & Piran (2000) argue for a ~ 0.1 \iG 
extragalactic field, corresponding approximately to equipartition between the magnetic, 
gravitational, and thermal energies of the Virgo supercluster for tangled fields. They also favor 
magnetic sheet topologies rather than two-dimensional random walk, and so a shorter distance 
results in the same deflection (~ 1 \iG would result in a lOn deflection for X, ~ 100 Mpc even at E 
> 6 X 10^^ eV). Observations will eventually reveal (an)isofropy, but modeling is called for to 
consfrain the possibilities for nearby GRBs. 

It is interesting that no similar Xiag-Xbreak-Qjet relationship has been established for the other 
major asfrophysical jet phenomenon, active galactic nuclei (G. Madejski, private communication). 
If a such relationship exists in AGN, its non-detection may be attributable to insufficient density 
of X-ray/y-ray observations in time and energy. For GRBs specfral lags can be determined on a 
relevant (and rapid) timescale that does not required scheduled, extended observations. Of 
course, to be more useful, this most fascinating phenomenon of GRBs will need to be observed 
simultaneously over a broader y-ray energy range. Future planned and proposed missions such 
as Swift, GLAST, and EXIST should provide very usefiil observations of GRB temporal 
behavior. Swift (Barthelmy 2001), with its improved sensitivity over BATSE and lower energy 
coverage, should give a larger yield of ulfra-low luminosity GRBs with immediate alerts, allowing 
searches for the corresponding supemovae. 

For many useful hints, suggestions, and encouragement I wish to thank: Jerry Bonnell, 
Jordan Camp, Joan Cenfrella, Thomas Cline, Dale Frail, Neil Gehrels, Kevin Hurley, John 
Krizmanic, Pawan Kumer, Greg Madejski, Miguel Morales, Robert Nemiroff, Jay Salmonson, 
and Jeff Scargle. Special thanks to David Band for a carefiil reading of the manuscript. 
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TABLE 1 

Quadrupole Moment in Supergalactic Coordinates for Long-Lag Bursts 



Tiag> Ngrbs Q Eq Na 



0.50 298 -0.025 0.019 L27 

LOO 155 -0.047 0.026 1.82 

1.50 96 -0.068 0.032 2.15 

2.00 72 -0.099 0.036 2.74 

2.50 54 -0.109 0.042 2.63 

3.00 41 -0.105 0.047 2.22 

3.50 27 -0.069 0.058 1.18 
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Figure captions 



Fig. 1 - Inferred afterglow beaming fraction vs. gamma-ray specfral lag for eight bursts with 
associated redshifts and BATSE data. The lag is corrected for time dilation, and approximately 
corrected for redshift of the spectral energy distribution (see section 2.4, Eq. [10]). Solid line, 
described by Eq. (la), was fitted only to the six bursts with 0.835 < z < 1.619; dashed line is toy 
model of Eq. (lb) for long-lag regime. Diamonds for beaming fractions derived by Frail et al. 
(2001); squares for those derived by Panaitescu & Kumar (2001). 

Fig. 2 - (a) BATSE peak flux (Fp, 50-300 keV) vs. CCF spectral lag (xiag , 25-50 keV to 50-300 
keV) for 1437 long-duration (T90 > 2 s) bursts. Dashed horizontal line indicates approximate 
BATSE trigger threshold. Solid lines divide the Fp-Xiag plane into three regions: bright short-lag 
bursts (region 1), dim short-lag bursts (region 2), and long-lag bursts (region 3). The point for 
GRB 980425 is circled, (b) The Xiag coordinate is magnified, showing only lags in the range -0.2 
< Tiag (s) < 1.0. Negative lags are completely consistent with larger errors at lower S/N. 

Fig. 3 - Ratio of numbers of bursts in regions 1+2 to region 3, Niong-iag / Nghort-iag > vs. peak flux 
binned in dyadic steps. Annotation shows that, while the sizes of both subsamples decrease near 
BATSE threshold, the increase in the relative frequency of long-lag bursts is significant. 

Fig. 4 - The integral peak flux-frequency disfributions (Log N[>Fp] - Log Fp) for the three 
subsample regions delineated in Figure 2. Notably, the long-lag bursts follow a -3/2 power-law 
over a factor of ~ 40 in peak flux. 

Fig. 5 - Sky disfribution in the Supergalactic coordinate system for very long-lag (xiag > 2 s) 
bursts (diamonds), suggestive of a tendency to concenfrate toward the Supergalactic Plane. Circle 
marks center of Virgo Cluster. Boxed diamond indicates position of GRB 980425; empty 
diamond is position of GRB 971208 (~ 20-30 s lag). Table 1 lists significances of quadrupole 
moment for different cuts in xiag. Q = -0.099 □ 0.036 for illusfrated distribution. 

Fig. 6 - Spectral lag distributions for peak flux ranges binned in dyadic steps. 

Fig. 7 - Similar to Figure 6, but with Xiag magnified to show the core region for short-lag bursts. 

2 1 

Overplotted histograms with gray fill are for brightest subsample (100 > Fp [photons cm s ] > 
8) S/N equalized to three lowest Fp ranges, demonsfrating that lower S/N level alone does not 
account for the broadening of the Xiag distribution at lower peak fiuxes. 
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Fig. 8 - Observed time-dilation factor vs. peak flux for T50 durations of those bursts in regions 
1+2, with fitted quadratic model. 

Fig. 9 - (a) The spectral lag distributions of Figure 6, overplotted (gray fill) with model lags 
including empirical errors, (b) Overplotted with the same model lags but without empirical 
errors. 

Fig. 10 - Similar to Figure 2a: scatter plot of Fp vs. model Xiag including empirical errors. 

Fig. 1 1 - Distribution of model lag with cosmological corrections for time dilation and spectral 
redshift for the two parameterizations of k in Eq. (10), ziim = 5 (solid) and 20 (dashed). 

Fig. 12 - Distribution of redshift, derived by varjdng z to match model peak flux to observed 
BATSE peak flux, assuming the two-branch relationship between xi^g and luminosity of Eq. (9); 
ziim = 5 (solid) and 20 (dashed). Abscissa, expressed as (1 + z), overemphasizes frequency of 
long-lag bursts, placing them in one bin. 

Fig. 13 - Distribution of proper distances associated with redshifts of Figure 12; ziim = 5 (solid) 
and 20 (dashed). Approximately 90 modeled burst sources lie at distances less than 100 Mpc. 

Fig. 14 - Scatter plot of peak flux versus redshift (1 + z). (a) ZHm = 5, (b) ziim = 100. The hard 
limit of Lcutoff is manifest in the upper right bound to the populated region. For ziim = 20 

sources actually pile up at this boundary, reflecting increasing redshift to accommodate short-lag, 

2 1 

high luminosity bursts (see text). Below Fp ~ 4 photons cm s , the concentration of low- 
luminosity, long-lag bursts appears near z = 0. 

Fig. 15 - Distribution of 24 redshifts determined to be associated with GRB afterglows or the 
host galaxy (see Barthelmy, GCN alerts: ht^://gcn.gsfc.nasa.gov/gcn/gcn3_archive.html). The 

mode and median continue to hold near z « 1, with ~ 2/3 of this sample within 0.4 < z < 2. 

2 1 

Compare to modeled redshifts in Figure 14b for bright bursts, Fp > 10 photons cm s . 
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Fig. 16 - Luminosity distribution in equal log-spaced intervals predicted by the two-branch lag- 
luminosity relation expressed in Eq. (9) for the modeled lags {zwm = 5) of the BATSE sample. 
The seven cascading dotted/solid histograms on right are for redshift-limited cuts: z < 30, 20, 10, 
5, 3, 2, and 1. In the regime z < 2 over the two decades 10"^'^ < L53 < 10^*^"^ the distribution is 

—1 8 

volume-limited, with fitted power-law slope of -0.82 (solid straight line), so dNyoi/dL L ' . 
For z < 1 within the range 10"^'^ < L53 < 10"^ the distribution is sensitivity limited, with power- 
law slope of +0.01 (dashed straight line), so dNsen/dL L~\ Sources with L53 < 10~^'^ 
(leftmost dotted histogram, z < 0.024, d < 100 Mpc) comprise the subsample with xiag < 2 s. 

Fig. 17 - Beaming fraction distribution for modeled lags of Figure 11 (ziim = 5) divided by fb, 
yielding distribution of all sources within BATSE reach, irrespective of beaming direction. This 
is still a sensitivity-limited disfribution, with the volume sampled decreasing as fb increases. 
Implied rate of highly luminous bursts (fb _ 3 x 10"^, L > 10"^^L53) is Ngrb _ 2 x 10^ yr~\ 

Fig. 18 -For profiled jet scenario: GRB sources within BATSE sensitivity reach irrespective of 
jet axis orientation with respect to observer. Luminosity dependence is L(0v) = 3.6 L53 x 
[0v/0o]^ with cutoff L > 10~^'^L53, for two values of the ininimvim jet cone radius, 0o = 1.50 
(solid histograms) and 3^ (dashed histograms), corresponding to 0.026 and 0.052 radians, 
respectively. In (a) y = 2.5, in (b) y = 2.0. Accounting for sample and BATSE exposure results 
in 0.66-0.85 x 10^ yr"^ (eo=1.5Q and 1.7-2.2 x 10^ yr"^ (©0=30 for the most highly luminous 
GRBs (L > 10-^^L53). 
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